Spin Liquid States on the Triangular and Kagome Lattices: A Projective Symmetry 

Group Analysis of Schwinger Boson States 
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A symmetry based analysis (Projective Symmetry Group) is used to study spin liquid phases on 
the triangular and Kagome lattices in the Schwinger boson framework. A maximum of eight distinct 
Z2 spin liquid states are found for each lattice, which preserve all symmetries. Out of these only 
a few have nonvanishing nearest neighbor amplitudes which are studied in greater detail. On the 
triangular lattice, only two such states are present - the first (zero-flux state) is the well known 
state introduced by Sachdev, which on condensation of spinons leads to the 120 degree ordered 
state. The other solution which we call the 7r-flux state has not previously been discussed. Spinon 
condensation leads to an ordering wavevector at the Brillouin zone edge centers, in contrast to the 
120 degree state. While the zero-flux state is more stable with just nearest-neighbor exchange, 
we find that the introduction of either next-neighbor antiferromagnetic exchange or four spin ring- 
exchange (of the sign obtained from a Hubbard model) tends to favor the 7r-flux state. On the 
Kagome lattice four solutions are obtained - two have been previously discussed by Sachdev, which 
on spinon condensation give rise to the q — and y/3 x \/3 spin ordered states. In addition we find 
two states with significantly larger values of the quantum parameter at which magnetic ordering 
occurs. For one of them this even exceeds unity, Kc « 2.0 in a nearest neighbor model, indicating 
that if stabilized, could remain spin disordered for physical values of the spin. This state is also 
stabilized by ring exchange interactions with signs as derived from the Hubbard model. 



I. INTRODUCTION 

Recent experiments on frustrated quantum magnets 
with unusual propertiesiiSi^ has revived interest in spin 
liquid phases. At the same time the theoretical under- 
standing of spin liquid phases in Z? > 1 has also ma- 
tured over the last few years - for example it is now es- 
tablished that such states are typically accompanied by 
emergent gauge fields in their deconfined phase. Sim- 
ple and tractable lattice Hamiltonians that realize these 
phases can also be constructed. Nevertheless, there is 
still a pressing need to understand the situations in which 
microscopic spin models might exhibit such phases, and 
the variety of such phases and their phenomenological 
properties. 

The triangular lattice antiferromagnet has been con- 
sidered a good place to look for such spin liquid states, 
its geometry is considered less conducive to magnetic or- 
der than for example the square lattice antiferromag- 
net, and more recently there has been evidence from 
numerical calculation s^i^i^i^ especially in the presence 
of anisotropic interactions and ring exchange, for un- 
conventional physics. Similarly, recent experiments on 
quantum magnets with a triangular lattice structure - 
Cs2CuCl4-\ K-(ET)2Cu(CN)32 and ^He films adsorbed 
on graphite^ show unusual properties that may be in- 
terpreted in terms of spin liquid phases or proximity to 
such a state. While the spin gapped Z2 spin liquid states 
considered in this paper are probably not appropriate 
to directly describe these experimental systems, except 



perhaps at the very lowest temperatures, they motivate 
a closer study of possible spin liquid states on the trian- 
gular lattice. The highly frustrated Kagome lattice has 
long been considered a good system for realizing exotic 
ground states. While there are fewer experimentally well 
studied examples of quantum spins on the Kagome lat- 
tice, conventional semiclassical analysis fails to produce a 
unique ground state when applied to this frustrated lat- 
tice. Hence, spin disordered states like spin liquids may 
be realized. 

In describing spin liquid states of quantum magnets, 
two distinct approaches have been found to be useful. 
The first involves writing the quantum spin (e.g. spin- 
1/2) as a fermion bilinear (fermionic spinons), and im- 
posing a constraint on the total number of fermions on 
each site. This fermionic approach is analytically pur- 
sued by writing down a mean field ansatz where the con- 
straint is imposed only on average to begin with. This 
mean field theory may also be viewed as a source of vari- 
ational wave- functions, which after projection yields spin 
wave functions. An important step in understanding the 
variety of spin liquid states possible in this representa- 
tion was taken by Wen and coUaboratorsSiiSiii who in- 
troduced the Projective Symmetry Group (PSG) classi- 
fication of spin liquids that respect the microscopic sym- 
metries of the system. This provides a guide to identify- 
ing the correct gauge group describing the spin liquid (if 
it survives projection) and also allows for a further clas- 
sification of different spin liquid phases with the same 
gauge group. These different spin liquid phases differ in 
subtle ways from one another, for example in terms of the 
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location of the wavevector of their lowest energy spin-1 
excitation. The advantage of such classification is that 
it allows for an exploration of possible spin liquid states 
independent of specific Hamiltonians that might realize 
them as ground states. In this way the symmetric spin 
liquid states within the fcrmionic spin representation on 
the square and triangular lattices were classifiediS*ii. 

A different way to describe spin liquids is via the 
Schwinger boson approach, where the spin operator is 
written as a boson bilinear (bosonic spinous) and was 
analyzed using a large-N approacbiSiiiiiii^ . An advan- 
tage of this approach is that it is able to readily access 
both spin liquid states as well as conventional magneti- 
cally ordered states, which arise from the condensation of 
the Schwinger bosons. For instance, in the case of the tri- 
angular latticei^, Sachdev found a spin liquid state which 
on decreasing the strength of quantum fluctuations to the 
size expected for s= 1/2, gave way to the 120 degree mag- 
netically ordered state, expected of classical spins on this 
lattice. However the possibility of other Schwinger bo- 
son spin liquid states on the triangular lattice have not 
been systematically explored. Similarly, on the Kagome 
lattice, SachdeAii^ found two solutions, which on spinon 
condensation gives rise to \/3 x \/3 or the q = magnet- 
ically ordered states. 

Mapping out the possible spin liquid states on these 
lattices is of interest for two reasons. First, when at- 
tempting to identify a candidate spin liquid state to 
match with experimental data or numerical simulations, 
one needs to specify precisely what the nature of the state 
is, and the quantum numbers of its excitations. There- 
fore while the simplest spin liquid may fail to match these 
details, other spin liquid states may give different predic- 
tions. Second, there may exist solutions stabilized with 
different interactions which remain spin liquids even for 
physical values of the spin. This would be particularly 
appealing as it would also provide direction to the search 
for spin liquid states. Even when the system is in the 
ordered phase, if it is proximate to a critical point where 
spin liquid physics holds, it has been arguedii^ that finite 
energy and finite temperature signatures of this will be 
unusual and best described in terms of the spin liquid 
variables. Therefore, a systematic investigation of the 
possible solutions of the Schwinger bosons states is re- 
quired on the lines of the PSG analysis of the fcrmionic 
spin liquid states. This will also shed light on the connec- 
tion between two approaches, for example, which states 
can be described in one but not the other approach. 

In this paper we will adapt the technique of Pro- 
jective Symmetry Groups, developed by Wen and 
collaboratorsSiiSiii in the context of fcrmionic mean field 
theories to study Schwinger boson mean field theories. 
We will focus on the triangular and Kagome lattices, in 
particular on spin liquid states with the Ising gauge group 
(Z2 gauge theories) which preserve the microscopic spin 
symmetries. Surprisingly, the number is not large (less 
than or equal to 8 for both lattices), which includes, in 
addition to the solution found by Sachdev, a set of fun- 



damentally different mean field wave functions, which 
give rise to Z2 spin liquids. In particular, if we make 
the further reasonable assumption of restricting to states 
with non-trivial nearest-neighbor bond amplitudes, there 
is only one other state on the triangular lattice, which we 
call the TT-flux state. Similarly, on the Kagome lattice, 
Z2 spin liquid states with nonvanishing nearest neighbor 
bond amplitudes, include, in addition to the two states 
found by Sachdev, two additional states, which we re- 
fer to as the [tt Hex,7r Rhom] and [0 Hex,7r Rhom] (which 
specifies the flux in the length-six hexagonal loop, and the 
length-eight rhombus). One of the problems in search- 
ing for spin liquid states has been the tendency of vari- 
ous spin Hamiltonians to yield magnetically order states, 
which is reflected in the small values of Kc, the critical 
quantum parameter at which spin ordering occurs in the 
Sachdev states: Kc ~ 0.3 on the triangular, and Kc ~ 0.5 
for the Kagome states in large-Ni^. (Note, k = 25" in 
the case of spin S but it is convenient to extend it to 
take on all real values). A feature of the new states that 
is their relative stability against magnetic order which is 
seen in the larger critical quantum parameter at which 
spin ordering occurs: Kc ~ 0.75 for the pi-flux triangular 
state, and Kc ~ 0.9 and Kc ~ 2.0 for the two new Kagome 
states respectively, in nearest neighbor models. 

While the distinct mean field solutions of the PSG 
analysis are typically local minima of the mean field en- 
ergy, an important question is whether they can be sta- 
bilized as global minima with appropriate interactions. 
This is discussed in detail for the case of the triangular 
lattice in this paper. While the zero-flux state is sta- 
bilized with just nearest-neighbor antiferromagnetic spin 
couplings, introducing next neighbor antiferromagnetic 
couplings or ring exchange (of the sign arising from the 
Hubbard model) both tend to favor the 7r-fiux state. We 
obtain the conditions on the microscopic spin interactions 
which favor this state and argue that this may be a fruit- 
ful parameter regime for numerical searches for spin liq- 
uids. In general we find that ring exchange on even length 
loops (with the sign as derived from the Hubbard model), 
favors states with tt fiux on length-4n loops, while it fa- 
vors states with zero flux in length- (4n + 2) loops (n is 
positive integer). This is in contrast to fcrmionic case 
where ring exchange disfavors fluxii for fermions in all 
loops. Addition or ring exchange on length-6 and length- 
8 loops can stabilize the [0 Hex,7r Rhom] state on the 
Kagome, which is good spin liquid candidate since it has 
Kc = 2.0. 

On spinon condensation, the triangular lattice 7r-fiux 
state naturally leads to magnetically ordered states with 
wavevectors at the midpoints of the Brillouin zone edges 
(in contrast to the zero-flux state which leads to the 120 
degree state with wavevector at the zone corners). This 
allows us to understand semi-classical (large spin) cal- 
culations in the presence of moderate-neighbor antifer- 
romagnetic couplings or " antiferromagnetic" , (with sign 
as obtained from the Hubbard model) ring exchange on 
the triangular lattice, where ordered states with the same 
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wavevector are found. It is tempting to connect the pos- 
sible ordered states on the triangular lattice as emerging 
from spin condensation out of the few Schwinger boson 
spin liquid states allowed by the PSG. New quantum 
transitions are expected on spinon condensation out of 
the spin liquid states obtained here, and will be the sub- 
ject of future study. 

Layout of the paper: In Section^we briefly review the 
Schwinger boson mean field theory. Section IIIII analyzes 
possible spin liquid states on the triangular lattice. It 
first reviews the Projective Symmetry Group classifica- 
tion of spin liquid states and the strong constraints that 
arise from relations between symmetry group elements. 
This is then applied to Schwinger boson states on the 
triangular lattice. A new state is found, the 7r-flux state, 
which is further analyzed - in particular spin configura- 
tions resulting from spinon condensation are described, 
and Hamiltonians stabilizing this mean field solution are 
obtained. The general effect of ring exchange interactions 
on Schwinger boson mean field states is discussed. In Sec- 
tion possible spin liquid states on the Kagome lattice 
are studied, and the properties of one of them, which 
is unusually stable against spin ordering, is described in 
more detail. The PSG analysis and other details are rele- 
gated to the appendices, which also contains analysis for 
other lattices of interest such as the anisotropic triangu- 
lar lattice. 



II. SCHWINGER BOSON MEAN FIELD 
THEORY 

There are a variety of ways of formulating the 
Schwinger boson mean field theory, for example, as a 
large-N approachi2*i^, or as an approximate variational 
approach. 

Here we will formulate a variational approach that will 
provide us with a unified way to study the effect of dif- 
ferent interactions. We write the spin Hamiltonian: 

if - Ji ^ S, • S, + J2 ^ S, • + • • • (1) 

in terms of Schwinger bosons: 



Si • Sj = -h\^(7„^„,bi„, ■ bj^ar.r'bjT' (2) 
with the constraint that at every site: 

T.^lb^.^n (3) 

cr 

where for a spin system with spin S, k = 2S. In the 
analysis below, it will be convenient to consider k to be a 
continuous parameter, taking on any non-negative value. 



We now consider a variational approach to finding the 
ground states and excitations of Motivated by the 
operator identity 

S. • S, =: Bj^B,, : -^i,, (4) 

where : : is normal ordering, and operators A and B are 
defined as 

(7 

Aij = ^'^ecrcr'btabja' (6) 
cr.cr' 

we consider a "Mean Field" Hamiltonian which is 
quadratic in terms of the Schwinger bosons, 

Hmf = J^J i~ A*^ A,, + B*jB,j +h.c.) 

i'j 

+ E Jij {A*jAij — B*jBij) 

where complex numbers Atj — ~Aji , Bij — B*^ are the 
parameters of the mean field ansatz. 

In the large-N Sp{N) theory, the mean field Hamil- 
tonian contains only the large-N generalization of the A 
term. However, since both A and B terms are consistent 
with global SU (2) symmetry (global spin rotation sym- 
metry) they are both included in the current theory. The 
introduction of both terms can be found in Gazzai^ and 
many other papersiS*^. 

This Hamiltonian is used to generate a variational 
wavefunction in terms of the variational parameters 
1 5* (Ay, Bij,fj,)). In order to obtain a spin wavefunc- 
tion, we need to project j^*) into the constrained Hilbert 
space where the total number of bosons at each site is 
exactly 25". Strictly speaking one must evaluate varia- 
tional energies after this projection step, using the spin 
Hamiltonian This generalization of Gutzwiller pro- 
jection to Schwinger boson has been studied by Chen and 
collaborators^iiSS. However, since this hard projection is 
not possible to implement analytically (it is even difficult 
to do numerically) we rely on a approximate strategy that 
forgoes implementing the constraint locally, but only on 
the average - i.e. we tune fi to ensure that: 

E(C^-) = « (8) 

(7 

we then evaluate the expectation value of the Hamilto- 
nian written out in terms of Schwinger bosons Q 
using the variational wavefunction. The resulting vari- 
ational energy is then minimized with respect to the 
variational parameters Aij , Bij . This yields the self- 
consistent equations, 

(iy)=Ay, (9) 
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When these self-consistent equations are satisfied, the 
variational(mean field) energy is simply obtained by uti- 
lizing the following identity, 

{S.-S,)^^{\B,,\'~\A,\') (10) 

The effect of projection on energetics will be a topic of 
future study. 



III. PSG FOR SCHWINGER BOSON STATES 
ON THE TRIANGULAR LATTICE 

We would like to classify Schwinger boson mean field 
states available to us on the triangular lattice in terms 
of the distinct spin liquid phases they can give rise to. 
This exercise allows us to readily obtain potentially in- 
teresting states relevant to quantum spin systems on the 
triangular lattice. Interactions that stabilize these states 
are discussed subsequently. 

This classification follows the classification of fcrmionic 
mean field states inSiiSiii. Here however the underlying 
local transformation is U{1) (not SU{2) as in the case 
of fermions). For mean field theories that survive pro- 
jection (or equivalently, survive fluctuations) to give rise 
to a spin liquid state, this procedure classifies the dis- 
tinct quantum phases of the system. We will require 
that the mean field theory refiects the underlying micro- 
scopic symmetries of the spin model. This leads to sym- 
metric spin liquid states. The symmetry transformations 
include the space group (lattice translations and point 
group symmetries) of the triangular lattice, spin rotation 
symmetry and time reversal symmetry. The novel ingre- 
dient here is that the mean field state may preserve these 
symmetries in some indirect way. In addition to the sym- 
metries of the original spin model, one also finds an extra 
global symmetry in all the mean field ansatz which is a 
subset of the local U{1) transformations mentioned pre- 
viously. This is called the Invariant Gauge Group (IGG) 
in the terminology of Ref?^. Since we are interested in 
mean field states that have exactly the same symmetry 
as the underlying spin model, not less or more, this is to 
be identified with the gauge group of the emergent gauge 
theory. In the examples below we will find an extra Z2 
symmetry in the mean field ansatz, hence the spin liquids 
that are obtained with this procedure are Z2 spin liquids, 
but with different internal structures depending on how 
the microscopic symmetries of the lattice model are real- 
ized. Thus, in contrast to conventional states which are 
distinguished by patterns of broken symmetry, spin liquid 
phases that are completely symmetric and even share the 
same gauge group can be further distinguished in terms 
of how the symmetries are realized. 

In Schwinger boson representation of spins, there is a 
local U{1) transformation of bosons: 

5r, ^ e'^^'^^h, (11) 



which leaves all physical observables unchanged. Under 
this transformation, our mean field ansatz (Aij, Bij) also 
transform: 

Ay -> e-"^'-'^-"''^^^A,j (12a) 
By -> e+^'t'^''^-"''^^^ B,j (12b) 

Two mean field ansatz that are related by such a transfor- 
mation, give rise to the same spin wave function after pro- 
jection. However, we will be interested in a class of trans- 
formations that leave the mean field ansatz invariant. 
Naively, one might expect that since we are interested in 
states that maintain all the microscopic symmetries, the 
mean field ansatz should be invariant under their opera- 
tion {e.g. lattice translations). While this is a sufficient 
condition for invariance, it is not a necessary one, due to 
the presence of the local U{1) transformations described 
above. A symmetry operation might return the ansatz to 
a U{1) transformed form which would suffice, since the 
same spin wavefunction will be obtained on projection. 
Hence symmetry transformations that leave the Ansatz 
invariant in general will contain the naive transformation 
combined with a local U{1) transformation. The set of 
all transformations that leave a mean ansatz invariant is 
called the Projective Symmetry Group (PSGf* . In princi- 
ple, we would like to associate each and every element of 
this group with a physical symmetry. This is because if 
the mean field ansatz is to be a faithful representation of 
the microscopic physics, it should have exactly as much 
symmetry as the original model. However, we always find 
that there are some elements of the PSG that are pure 
local transformation of the kind Hll|) . The set of such el- 
ements also forms a group (a subgroup of the PSG) and 
is called the invariant gauge group (IGG). These cannot 
be the result of a physical symmetry. It is therefore nat- 
ural to associate these elements with the emergent gauge 
group that describes the spin liquid phase obtained (if 
the mean field state survives projection). Therefore the 
first step is to identify the IGG of a mean field ansatz. 

The Invariant Gauge Group for Schwinger Boson 
Mean Field Theories: A general Schwinger boson mean 
field Hamiltonian with explicit global SU{2) symmetry 
(global spin rotation symmetry) must be of the form of 
equation |7J). 

It is clear that if Aij and Bij are both nonzero the IGG 
must be a Z2 group. The only two elements of IGG are 
identity operation 1 and the IGG generator —1 : b^a 
—bra- This can be seen by considering only one bond. 

Note, that if Bij are nonzero while all Aij vanish, the 
IGG will be a U{1) group: b^a e^^bro-, where cj) is 
a site-independent constant. However it is unlikely for 
this ansatz to describe an antifcrromagnet according to 
(|10|l . Finally we should consider the case with only non- 
vanishing Aij. On a frustrated lattice the IGG will still 
be the above Z2 group. However, if the lattice is bipar- 
tite, the IGG will be a U{1) group: b^a e'^^'^bra, where 
we apply opposite signs on the two sublattices. This is 
the case of simple square lattice. 
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We will only consider PSGs with Z2 IGG in the follow- 
ing discussion. Hence we are implicitly restricting our- 
selves to Zi spin liquid states, which are the natural spin 
liquid states on frustrated lattices within the Schwinger 
boson formalism. 



A. Algebraic Constraints on the PSG 

We consider mean field Hamiltonian that preserves all 
of the physical symmetries in the PSG sense (symmetric 
spin liquid states). They are spin rotation symmetry, lat- 
tice space group symmetries and time reversal symmetry. 
The spin rotation symmetry is already implemented by 
considering mean field ansatz of the form , which is ex- 
plicitly invariant under global SU{2) spin rotations. We 
will consider time reversal symmetry at the end of our 
derivation. The operations that we will now pay special 
attention to are the space group symmetries (translations 
and point group operations for the triangular lattice). As 
discussed before, these can be implemented via combin- 
ing the naive transformation with a local (gauge) U{1) 
transformation. One can ask the question - is it pos- 
sible to have mean field ansatz with any choice of the 
gauge transformations? It turns out that there are alge- 
braic relations among the symmetry group elements that 
strongly constrain the possible choices of gauge transfor- 
mations. Thus, to obtain all possible PSGs we should 
first check the algebraic structure of PSGs without refer- 
ence to a specific mean field ansatz. The possible PSGs 
allowed by the algebra of the space group are defined as 
the algebraic PSGs related to the space group. 

For the isotropic triangular lattice, the space group is 
generated by two translations Ti and T2, one reflection a 
in a bond and the 60 degree rotation R:^ about a lattice 
site. 

We use the following oblique coordinate system. In 
this system a site index r has two integer components 
r — (ri,r2). Then the four generators are given by the 




FIG. 1: Oblique coordinates system of triangular lattice. 



such that the mean field ansatz is invariant under Gti?!, 

GT2T2, GaU, and Gr^Ril, respectively. We can repre- 

3 3 

sent the four gauge operations by their phase. 



G 



X 



J(t>x{r)i 



(14) 



where the X are Ti, T2, cr, and R:^, respectively. 

The PSG is generated by combining the generators of 
the IGG: —1, and the above four compound operators 
GxX. 

The structure of the space group imposes algebraic 
constraints on the Gx- For instance, the combination of 
translations shown in FIG. |21 should be equivalent to the 
identity, i.e. Ti~^T2TiT2~^ = T. Therefore we require 
that the implementation of these symmetries in the PSG: 
(GtJ^xY^Gt2T2GtJ^y{Gt2T2)~^ must be the equivalent 
of an identity operation, which means it is an element of 
the IGG, either 1 or —1. This string can be rewritten 
as Ti-^Grr^Ti ■ T^-^Gt^Ti ■ {T^-^T2)GtATi~^T2)-^ ■ 
(Gts)"^- Using the fact that for a space group operation 
Y , the gauge transformation {Y)^^ ■ Gx ■ Y acting on 
site r will just give a phase 0;is:[y(r)] (where Y{r) is the 
image of r under the space group operation Y), we end 
up with the equation 



=pi7r 



% [Ti (r)] + (/.T, [Ti (r)] + c/^t, [T2-'T, (r)] ~ (/-t, (r) 

(15) 



where pi € {0, 1} is independent of the site index r, and 
arises from the fact that the identity operation can be 
any one of the two elements of the IGG. 
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FIG. 2: An example of a relation between space group el- 
ements Ti~^T2TiT2~^ ~ ^, that leads to an algebraic con- 
straint on the PSG. 

This kind of algebraic constraints strictly restrict the 
possible structures of PSG. It turns out (proof is in Ap- 
pendix that there is only a finite set of such con- 
straints, which if satisfied guarantees that all relations 
between space group elements are satisfied. These rela- 
tions, beginning with the one described above, are: 



following formulae 






T2T1 


- T1T2 


(16a) 






TiRjL 




(16b) 


Ti ■■ {ri, ra) - 


^ (ri + 1, ra) 


(13a) 


T2R^ 


= R^TiT2 


(16c) 


T2 ■■ {ri, ra) - 


~> (ri, r2 + 1) 


(13b) 


Tin 


= (JT2 


(16d) 


o- : in, r2) - 


{r2, ri) 


(13c) 


Tin 


= (TT2 


(16e) 


Rf ■■ in, r2) - 


(ri - r2, ri) 


(13d) 


3 


= 


(16f) 


For symmetric spin liquid states there must be four 


a-' 


= a 


(16g) 


gauge group transformations G^i, Gjvj, G„, 


and Gj}„ , 
7 


RlLG 

3 


= 


(16h) 
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These relations place severe constraints on the allowed 
PSGs. The most general solution that satisfies them is: 

</>Ti(ri, ra) = (17a) 
</>T2(?'i, '"2) = PiTTi (17b) 
0(t(»'i, '"2) = _P27r/2 +pi7rrir2 (17c) 

2(/>i?,i (?■!, ?'2) = PsTT +Pi7rr2(r2 - 1 + 2ri) (17d) 

3 

where pi, p2, P3 are either or 1. Thus there are 
at most eight distinct Schwinger boson symmetric spin 
liquid states on isotropic triangular lattice. If we put 
more conditions on the mean field ansatz {e.g. that 
the nearest-neighbor Aij are nonzero), then the number 
of possible symmetric spin liquid states is reduced from 
eight. 

A detailed derivation of the above formulae is given in 
Appendix ^ Also we include the solution of algebraic 
PSGs on the anisotropic triangular lattice. 



B. From PSGs to Mean Field Hamiltonians : 
Nearest-neighbor Models 



real and nonzero in general. Explicit expression for this 
ansatz are: 

^(''l^''2)^(''l + l,''2) ^(ri ,^2)— ►(ri ,r2 + l) 
= ^^(ri,r2)^(ri + l,r2 + l) = ^li 

-S(ri,r2)^(i'i+l,r2) ^ ^ (ri ,r2)^{ri ,r2 + l) 
= -^(ri,r2)-^(ri+l,r2 + l) = Bl 

The PSG predicts that next-nearest-neighbor A must be 
zero. Since Aij = ^Aji is real and anti-symmetric, it 
is natural to represent it by oriented bonds. FIG. His 
a pictorial representation of the ansatz for the zero-flux 
state. 




If we assume that nearest-neighbor amplitudes Aij are 
nonzero, (which is natural given that nearest neighbor 
interactions in physical models tend to be dominant and 
usually antiferromagnetic) , there arc more constraints on 
the possible PSG structures. 

This can be seen as follows. Since a maps bond 
(0,0) ^ (1, 1) to itself, we must have (/)„(0, 0)-!-^^ (1, 1) = 
mod 27r. This imposes the constraint p2 ~ pi mod 2. 
Also, bonds (0,0) ^ (-1,0) and (0,0) ^ (1,0) are 
related by a 180 degree rotation, (R^)^, and by anti- 
symmetry and translation Ti we have A(o,o)^(-i.o) — 
0)^(0.0) = —A(o, 0)^(1,0)- This leads to another 
constraint 

</)fl^ (-1, 0) + (0, 1) + <^fl^ (1, 1) + 3^R^ (0, 0) 

3 333 

=7r mod 2n 

which fixes p^ = 1 — pi mod 2. 

Thus we have only two non-equivalent PSGs, corre- 
sponding to pi = 0, which we call the zero-flux state or 
Pi = I which we call the tt-Aux state. 

1. The zero- flux state: 

This state is specified in terms of the integers pi = 
0, p2 = 0, P3 = 1, or equivalently in terms of the phase 
factors involved in implementing the space group opera- 
tions: 

0Ti(fi, r2) ^ (f>T2{ri, r2) ^ ipain, r2) = 

(l)R^iri, rs) = it/2 

We can now ask what mean field ansatz would be char- 
acterized by such a PSG. The mean field ansatz is speci- 
fied by the amplitudes Aij and on the various bonds. 
For this PSG, Bij and Aij on nearest-neighbor are both 



FIG. 3: Ansatz of the zero-flux state. An arrow from site i 
to site j means Aij > 0. All Aij have the same magnitude. 
All Bij are real and of uniform magnitude (not shown). 

This ansatz has explicit translational invariance, and 
has been studied by Sachdev using the large-N method 
2. TT-flux state: 

This state is specified in terms of the integers pi = 
1, p2 = 1, J53 = 0, or equivalently in terms of the phase 
factors involved in implementing the space group opera- 
tions: 

0Ti(fl, f2) = 
4>T2{tI, f2) = TITl 

't'c^iri, r2) = (7r/2) + 7rrir2 

{ti, r2) = 7rrir2 + (7r/2)r2(r2 - 1) 

3 

The mean field ansatz that realizes this PSG is as fol- 
lows. While the nearest-neighbor bond Aij is real and 
nonzero, the nearest-neighbor bond Bij must be zero. 
Expression of the ansatz is: 

(— I)'"^^(ri,r2)^(ri+l,r2) — ~^(ri ,r2)^(ri ,r2-|-l) 
= — (— l)''^A(rj^r2)^(ri + l,r2 + l) =^1, 

-^('■l.r2)^(i-l + l,r2) — B{ri,r2)^{rx,r2 + 1) 
— -^(n,r2)^('-l + l,-r2 + l) = 

The PSG also predicts that the next-nearest-neighbor B 
must be zero. Note, for this state the translation symme- 
try is not explicit in the Mean Field Hamiltonian, hence 
the unit cell for the Schwinger bosons is doubled. This tt- 
flux ansatz for nearest-neighbor model is shown in FIG. 01 
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FIG. 4: Ansatz of the 7r-flux state. All nearest-neighbor Aij 
are real and of the same magnitude, and are positive in the 
directions shown. All nearest-neighbor Bij are zero. The unit 
cell for the Schwinger bosons is doubled due to the flux, and 
u and V are the two sites in this unit cell. 

These two states can also be distinguished 
by the gauge-invariant phase of the product 
Aij{~A*f.)Aki{~Ai^), where ij,kj form a rhom- 
bus with unit length sides. This quantity was introduced 
by Tchernyshyov, et. alm^ as the flux in the rhombus in 
bosonic large-N theory. Defining: 

|Ai|V* = A,,(-A*,)A,,(-A*) (18) 

where \Ai \ is the uniform magnitude of Aij, we have that 
for the zero-flux state, the flux $ = for all rhombi; while 
for the TT-flux state $ = tt for all rhombi. Therefore these 
two states are clearly not gauge equivalent mean field 
states. 

Finally we can consider time reversal symmetry T. 
Because the two ansatz are both real, they directly re- 
spect T-symmetry, since time reversal transformation 
will change the ansatz to their complex conjugate. It 
is then an interesting question to ask what kind of PSG 
can support T-breaking ansatz. It turns out that one 
must also break lattice reflection symmetry to obtain a 
time reversal breaking state. In Appendix. iBlwe list the 
solutions of algebraic PSGs for anisotropic triangular lat- 
tice and the realizations in nearest-neighbor model which 
do support a T-breaking ansatz. 

IV. ANALYSIS OF THE MEAN FIELD 
THEORIES 

We now study the mean field theories arising from 
the ansatz described previously. Given the PSG clas- 
sification, it follows that these different mean field solu- 
tions will be local minima of the mean field energy. In 
order to pick which of these is favored with a particu- 
lar Hamiltonian, one needs to compare energies between 
these states, which is left to the next section. Here we 
content ourselves with describing the properties of each 
of these mean field ansatz. The discussion closely par- 
allel Sachdev's large-N solution of triangular lattice and 
Kagome latticei^. We also take the quantum parameter 
K as a continuous number, although for Schwinger bosons 
derived from SU (2) spins of size S, k = 2S. For small k, 



the spinon dispersion will be gapped and we have a spin 
liquid mean field ground state. When k goes beyond cer- 
tain critical value Kc spinon dispersion becomes gapless, 
bosonic spinous condense, and magnetic long range order 
develops. For simplicity we present only the solution to 
the nearest-neighbor model. 



A. Zero-flux state 

The mean field theory of this state is almost identi- 
cal to Sachdev's large-N theory for triangular latticei^ 
and the phases obtained are continuously connected to 
the phases described in that work. The only difference is 
that we include the parameter Bij, which is not present 
in the large-N treatment, and therefore our spinon dis- 
persion and critical quantum parameter Kc differ slightly. 
This is reviewed briefly here before turning to a similar 
analysis of the tt-Aux state. 

The Brillouin zone for the triangular lattice and our 
choice of coordinate systems in k-space in shown in 
FIG. [S| In this oblique coordinates system, k • r = 
kiri + k2r2- For covenience we define = —ki — k2. 




FIG. 5: Triangular Lattice Brillouin Zone (BZ) and k- 
space coordinate system. The large hexagon is the BZ of 
the original lattice and the zero-flux ansatz. The large rect- 
angle is the BZ of the 7r-fiux ansatz. Black triangles at 
K-points: ±(27r/3, 27r/3) are the minima of the spinon dis- 
persion of the zero-flux state, and are also the wavevec- 
tors of magnetic ordering of the zero-flux state (120 degree 
state). Black diamonds at ±(7r/2, — 7r/2) are the minima of 
the spinon dispersion of the 7r-fiux state. Black hexagons at 
M-points: ±(7r, — tt), ±(— 7r,0), ±(7r,0) are the magnetic or- 
dering wavevectors of the 7r-flux state. 

After Fourier transformation, the nearest-neighbor 
mean field Hamiltonian (O is 

hmf= ;^^'(k)ti?(k)*(k) 

k 

+ N.lfi + fiK - MiiBf ^ Al)] 
where Ng is the number of sites and the vector spinon 
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field \E'(k) and the 2-by-2 matrix D{k) are: 

- (Hi) (M) 

n(\r\ - f JiBiRe{^k) - n JiAiIni(^k) A 
^^•"^ - zJiAiIm(Ck) JiSiRe(a)-My ^ ^ 

where Ck = e"=i + e^''^ + e-''-'''+''^\ 
After a Bogoliubov transformation 

where Mk e >S'C/(1, 1), is chosen to diagonahze the mean 
field Hamiltonian: 

Hmf = ^w(k)(7^^7kT +7lki7-ki + 1) 

k 

+N,[^i + ^iK~ 3Ji{Bl - Aj)] 
where the spinon dispersion a;(k) is: 



kya 




-4-2 2 



FIG. 6: The lower edge of the two-spinon spectrum for the 
zero-flux state at k = 0.3. Axes are in dimensionless unit 
kx,ya where a is lattice constant. Darker regions has lower 
energy, the 120 degree magnetic order arises from magnon 
condensation at the zone corners. See also FIG.|S] 



u;(k) = J [JiBiRe(ek) - /i]' - [JiAilm{(k)f (22) 



B. TT-flux state 



and the self-consistency equations are: 

6J1A1 = - / ^^^^^d^/c (23) 

(24) 
(25) 



Jbz dAi 

6J,B, = + / ^d^fc 
Jbz dBi 



1 + K 



9w(k) 2 



BZ 



where the integral is over the Brillouin Zone, and the 
integration measure d^k = dfcidfc2/(47r^). 
The mean field energy per bond is 



S/bond = Ji{Bf - Al) 



(26) 



Minima of the spinon dispersion are at (fci, ^2) = 
±(27r/3, 27r/3), which are the corners of the Brillouin 
zone. The two-spinon spectrum also has minima at the 
Brillouin zone corners and at the zone center. As in 
Refiii, we display in FIG. El a contour plot of the min- 
imum energy required to create a two spinon excitation 
at a given crystal momentum (the lower edge of the two- 
spinon spectrum) to show this feature. 

By solving mina;(k) = together with self-consistent 
equations, we determined the critical quantum parame- 
ter Kc ~ 0.42. We can estimate the sublattice magnetiza- 
tion for spin- 1/2 Heisenberg model to be 58% of classical 
spin, which probably overestimates the order since it is 
larger than that obtained from spin-wave theories which 
predicts this to be about 48%2i. 



From FIG. 01 it is clear that there are two sites in a 
unit cell, distinguished by the parity of r2- We label them 
by u and v. 

After Fourier transformation, the mean field Hamilto- 
nian O is 



H 



MF 



k 



'J'(k)^D(k)*(k) 



(27) 



Here sums over {Ns/2) k-points in spinon Brillouin 
zone, which is one half of the hexagonal Brillouin zone of 
original spin model (see FIG.[S1). 
The vector spinon field ^'(k) is: 



vE'(k) 



/ 6ukT \ 

^u-ki 

V^Lkx/ 



(28) 



and the 4-by-4 matrix D(k) has the following block form 



M-1 -^P0^)\ 

ip{ky /i • 1 J 



(29) 



where 1 is the 2-by-2 identity matrix, and the 2-by-2 
Hermitian matrix P(k) can be written in terms of the 
Pauli matrices a^'V^^: 

P(k) = JiAi [sin(/ci)CT^ - sm{k2)<T'' - cos{k3)ay] (30) 
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This Hamiltonian can be diagonalized by a Bogoliubov 
transformation using a SU(2,2) matrix after which we 
have: 

Hmf = ^ ^(k)(7lkT7akT +7l-ki7a-ki + 1) 

k; a—u,v 

+Ns[n + HK + 3JiAl] (31) 

where 7 are transformed boson operators and a — u^v is 
a sublattice index. 

The spinon dispersion is therefore four-fold degenerate 
with the form: 

Uj{k) = .y//^2 _ J2^2 [sin2(fci) + sin2(fc2) + cos2(fc3)] 

Self-consistent equations comes from minimization of 
mean field energy (constant terms after Bogoliubov 
transformation) . 



6J1A1 



1 



BZ 



aj(k)^2 



BZ 



(32) 
(33) 



where the integral is over the Brillouin Zone, and d^fc = 
dfcidfc2/(47r2). Recall, the amplitudes Bij on nearest 
neighbor bonds are forbidden in this state by the PSG. 

On solving the self-consistent equations, the mean field 
energy per bond is obtained as: 

£;/bond = -JiAj (34) 

and the minima of spinon dispersion are at (A;i,fc2) — 
±(7r/2,— 7r/2). Therefore the minima of the two-spinon 
spectrum are at the Brillouin zone edge-centers and cen- 
ter, as illustrated in FIG. [7| The distinction from the 
analogous diagram for the zero-flux state in FIG. |Blindi- 
cates that the 7r-flux state is really a distinct state from 
zero-flux state. 

The critical quantum parameter for 7r-flux state is 
found to be Kc ~ 0.75. This large Kc makes the 7r-flux 
state a promising candidate of spin liquid especially given 
that in the zero-flux case the mean fleld theory seemed 
to overestimate magnetic order. Even in the event that 
magnetic order is present, it is expected to be weak, and 
the presence of a proximate spin liquid should be appar- 
ent at finite temperatures and energies as discussed iiJiS. 
However, for the nearest-neighbor antiferromagnet, the 
zero-flux state always has lower mean field energy, which 
is consistent with the general argument of flux-expulsion 
by Tchernyshyov, et. al.^^ . Later we will see that tt- 
flux state can be stabilized as the global minimum of the 
mean fleld theory, if other terms such as next neighbor 
couplings or ring exchange are present in the Hamilto- 
nian. 



C. Spin Configurations from Spinon Condensates 

When quantum parameter n goes beyond its critical 
value the spinon dispersion becomes zero at several k- 




- kxa 



-4 



-2 



FIG. 7: The lower edge of the two-spinon spectrum of the 
TT-flux state with k = 0.3. Axes are in dimensionless unit 
kx^ya where a is lattice constant. Darker regions have lower 
energy. Magnetic order arises from magnon condensation at 
the centers of the zone edges. See also FIG.|S] 



points. Bosons condense at those k-points, namely, the 
spinon field ^'(k) receives expectation value in ground 
state. We briefiy analyze the structure of condensates 
and the corresponding spin configurations for the zero- 
fiux state, where we follow referencei^, and then apply 
the same analysis to the 7r-flux state on the triangular 
lattice. 

The zero-flux state, has spinon minima at kc = 
(27r/3, 27r/3). At this point, D{kc) of equation ((^OJ) has 
an eigenvector with zero eigenvalue 

(^ is transposition). When k just exceeds Kc, a conden- 
sate at kc is obtained (^(kc)) = ((6k^t) (^Ik^i)) ~ ^i^*!. 
Similarly, Z?(— kc) has one eigenvector associated with 
zero eigenvalue 

^2 = {-h if 

Condensate at — kc is (^'(— kc)) ~ C2'^2- Here Ci_2 are 
two complex numbers. 

The condensate on lattice site r is then given by 

ike' 



ZCi, 
C*2, 



-IC2 



and the ordered magnetic moment is S(r) = [XjT^x^ ax. 

This can be readily shown to give rise to the 120 degree 
classical Neel state, in which the classical spin is 



S(r) = ni cos(Q • r) + n2 sin(Q • r) 



(35) 



where r\.\^2 are orthogonal to each other and have the 
same length as the classical spin 

= = (36) 
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Q = ±(27r/3, 27r/3) is the wavevector of Brillouin zone 
corner, represented as small triangles in FIG.|S1 

For TT-flux state, the critical spinon wavevector occurs 
at kc — (7r/2, — 7r/2). At this point,, -D(kc) of equation 
has two eigenvectors associated with zero eigenvalue 



*i = -73' 0' 1 



*9 



\V3' Vl' ^' 



The condensate at kc can in general be parameterized 
as (*(kc)) = ci*i +C2*2- Similarly, -D(-kc) has two 
eigenvectors associated with zero eigenvalue 



*3 = 73' 0' 1 



*4 = (-73' -w' 1' o) 

and the condensate at — kj, is parameterized as: 
kc)) = €3^*3 + C4^'4. Here ci^2,3,4 are four complex 
numbers. 

Returning to real space, the spinon condensate on sub- 
lattice u is 



{bur^ 

{burl 

W^i + 73-^' 



V3 



C4 



while the spinon condensate on sublatti 



ice V IS 



_ i 
V3 



Cl 



1+i 

V3 



C2, 



V3 



C3 



V3 



C4 



^3: 



The ordered magnetic moment on sublattice p is 
(1/2) (xp)^ axp for p — u,v respectively and is present 
at wavevectors at the mid points of the Brillouin Zone 
edges shown via black hexagons in FIG. El 

Under the constraint that the magnitude of condensate 
is uniform (which is required to have uniform magnitude 
for the ordered moments ), we find a continuum of four- 
sublattices ordered states, which is consistent with previ- 
ous classical energy analysis on triangular lattice magnets 
with interactions that favor order at these wavevectoraS^. 
The classical spin in the four-sublattices states is 



S(r) = n,{-lY' + n2{-lp + n3{-lY'+^^- 



(37) 



where ni_2,3 are three vectors orthogonal to each other, 
and 

nl + nl + nl = S^. (38) 

Generally speaking this configuration is non-coplanar 
if I Hi ^2,3 1 are all nonzero. However, previous analysis con- 
firmed the existence of " order from disorder" phenomena. 



namely, quantum or thermal fluctuation will lift the ac- 
cidental degeneracy and favor a coUinear state in large 
S models with next nearest neighbor antiferromagnetic 
exchange^SiSi, which corresponds to the case that only 
one of ni_2,3 is nonzero. 

It has been argued that ring exchange will favor a non- 
coplanar stated. The classical ground state will then be 
a tetrahedral configuration with equal angle between any 
two of the four spins. In the above expression, this state 
is realized when ni_2,3 have same magnitude. 

However, on studying the ordered state with weak 
spinon condensation, the constraint that the magnitude 
of the spinon condensate is uniform enforces the addi- 
tional condition on the ni 2,3, that one of the following 
three equations must be satisfied 



|ni| = |n2| + IngI 
or |n2| = IngI + |ni| 
or |n3| = |ni| + |n2| 



(39a) 
(39b) 
(39c) 



Thus the classical degeneracy is not entirely lifted at 
this level. Also, neither the coUinear state nor the tetra- 
hedral state can be obtained from weak spinon conden- 
sates from this spin liquid state - although the wavevec- 
tors of all these ordered states are identical. This implies 
that if there is a continuous phase transition out of the 
TT-flux spin liquid phase, the proximate spin ordered state 
is not the coUinear or tetrahedral state, but these are re- 
alized via further phase transitions on increasing the size 
of the spin. We are still seeking a simple explanation for 
these results. 



V. HAMILTONIANS STABILIZING THE 
TT-FLUX STATE 

The TT-flux spin liquid state differs in several respects 
from the zero-flux state; also, its stability against mag- 
netic ordering up to a relatively large quantum param- 
eter Kc = 0.75 (in a nearest neighbor model) suggests 
its physics may be important for understanding physi- 
cal spin- 1/2 models. Hence it is of interest to ask what 
interactions might stabilize such a Tr-flux state. 

A clue is provided by the fact that the lowest en- 
ergy triplet excitations of such a state are located at 
the midpoints of the Brillouin zone edges. This is pre- 
cisely the ordering wave-vector for classical spins on 
the triangular lattice in a model with next-nearest- 
neighbour ( J2) antiferromagnetic interactions in addition 
to nearest- neighbour antiferromagnetic interactions (Ji) 
in the range 0.125Ji < J2 < Ji^^. It is conceivable that 
on reducing the size of the spin in this model from in- 
finity (the classical limit) to S=l/2, the system enters a 
spin liquid state described by the tt-Sux state, although 
strictly the quantum parameter for this case is k = 1. 

A second possibility arises from ring exchange - it 
is well known that there are several realization of 
S=l/2 quantum antiferromagnets on the triangular lat- 
tice where ring exchange plays an important role. These 
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include monolayers of "^He adsorbed on graphite^, and the 
organic quantum antiferromagnet /v-(ET)2 Cu2(CN);j2>. 
Also, exact diagonalization studies of the triangular lat- 
tice S=l/2 magnet with ring exchange has uncovered 
unusual properties^. In a recent variational study with 
fcrmionic mean field states by Motrunich^ it was found 
that ring exchange of the sign that arises naturally with 
electrons, favors a state with zero flux. Here we will 
investigate the effect of ring exchange on the Schwinger 
boson states - in particular for which sign of the exchange 
one might favor the 7r-fiux state. We find that including 
four spin ring exchange with the sign as obtained from 
the Hubbard model stabilizes the tt-Aux state. 

Here we give some physical arguments as to why next- 
neighbor antiferromagnetic interactions and the ring 
term will stabilize 7r-flux state over the zero-flux state. 



A. Effect of Next Neighbour Interactions J2: 

For zero-flux state, nearest-neighbor amplitudes Bi 
and Ai, and next-nearest-neighbor amplitudes B2 are 
nonzero, while the next-nearest-neighbor amplitude A2 
is zero. Note that B2 is nonzero even if next-nearest- 
neighbor coupling J2 is zero. 



zero— fiux 



bond • Ji 



B 



zero— fiux] \ ^ / ^[zcro— flux] ^ ^ 



. J2 / T-> [zero— fiux] ^ ^ 

In contrast, for the 7r-flux state, Ai and A2 are nonzero, 
while Bi and B2 are zero. Note, A2 is nonzero even if J2 
is zero. 



■fluxl 



bond • Ji 



. [tt- fiux] 



Ji 



A 



[tt— flux] 



It is readily shown with the self consistent values of 
the mean field amphtudes that for J2 = 0, < 
^[tt-Aux] ^ j ^Yie zero-flux state is preferred at the mean 
fleld level. However, if we increase J2 from zero, we do 
not expect the ansatz to change rapidly, thus the en- 
ergy for zero-flux state will have a positive slope with 
respect to J2, while the 7r-flux state has a negative slope. 
We may then expect a first-order transition between the 
two states. This is consistent with previous analysis 
of Ji — J2 model in the opposite limit of semiclassical 
spins, where a transition between the 120 degree state 
(the magnetically ordered analog of the zero-flux state 
and the colinear state, the magnetically ordered analog 
of the TT-flux state is found on increasing J2^^'^*- From 
the classical energy analysis^^, we know that in large- S" 
limit a = 1/8 is the critical point between different or- 
dered states with ordering wavevectors at the BZ corners 
vs. BZ edge midpoints. Thus one might expect that a 
moderate next-nearest-neighbor term (1/8 < a < 1) will 
stabilize the tt-Aux state over the zero-flux state in the 
quantum limit, since these states have short range order 



(two spinon minima) at precisely these wavevectors. We 
now provide estimates for the phase boundary (the crit- 
ical ratio a = J2/J1) between the zero-flux and 7r-flux 
spin liquid phases as a function of the quantum parame- 
ter K. 

Small K analysis: The competition between different 
spin liquid states can be studied analytically in the limit 
of small-K^^. In this case the magnitude of the ansatz 
Aij and Bij are small. Then, we can develop a series 
expansion of the self-consistent equations in k, and find 
solutions. 

First, we use this expansion to look at the J1-J2 model, 
with a = J2/J1 < 1 and sufficiently small quantum pa- 
rameter K. 

The energy difference between the zero-flux and the 
TT-flux state in this limit is: 



3-fli; 



E^ 



r-flu 



bond • Ji 



(1 



a 



'Hi' 


2 


2a 


' K' 


.6. 




1 -a 


.6. 



— + 23a -t- 72a^ 



87 

4 





' K' 


) 


.6. 



(40) 



From the above expressions for mean fleld energy dif- 
ference, we can find the transition point between the zero- 
flux and TT-flux states to lowest order in k. 



1 



(41) 



In the limit of vanishing k the zero-flux state is preferred 
until the next-neighbor bond is stronger than the near- 
est neighbor bond. However, this critical ratio decreases 
rapidly with increasing k. 

In the small k limit one can ignore the spin carrying 
excitations. Then, going beyond the mean field theory 
the relevant degrees of freedom are the gauge excitations, 
which in this case correspond to those of an Ising gauge 
theory. The distinction between zero and TT-flux states in 
this limit we believe has to do with the the sign of plaque- 
tte energy term that penalizes Ising vortex excitations in 
the case of the zero-flux state, but prefers a ground state 
with a uniform background "flux" of such vortices in the 
case of the TT-flux state. 

Numerical Solution: The mean fleld equations can be 
numerically solved to obtain the competition between the 
zero-flux and Tr-flux states and states with magnetic or- 
der. This is shown in FIG.|H1 The critical value k}c~^^^^ 
at which spinon condensation occurs for the Tr-flux state 
is shown by the dashed and solid red line. If the Tr-flux 
state is stabilized, then above this line spinon condensa- 
tion leads to magnetic order at the B.Z. edge centers (M- 
points). There are a number of distinct ordering patters 
consistent with this ordering wavevector - which makes 
an analysis of the ordered phase a delicate one; and we do 
not attempt it here. Therefore, the phase diagram shown 
in FIG.IHlis strictly speaking only to be trusted below this 
^[tt flux] jjj^g^j.g^ line). However, certain extrapolations 
above this line can be made with some degree of confi- 
dence. 
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The first-order phase boundary between the tt-Aux and 
zero-flux spin liquid states is shown by the black line in 
FIG. IHl At small k it approaches a = 1 via the analytic 
form derived previously 1411 It intersects the spinon con- 
densation line for the 7r-flux state at a = 0.46, k — 0.47. 

The zero-flux state is stabilized at smaller values 
of a, and sufficiently small quantum parameter k < 
^[zero flux]^ where the critical quantum parameter for 
spinon condensation in this state is shown by the blue 
line. Spinon condensation is expected to lead to the 120 
degree state on crossing this line-- . Again, this is most 
reliably established if the ordered states with wavevec- 
tors at the M-points (M-states) are not in the picture, 
which is the case below the dashed red line in FIG. |H1 
An interesting aspect of the phase diagram is the large 
values of ^1.^'='°"^"'^] ~ i_g i^j-^^t arise, if the M-states are 
neglected. This is believed to be a result of frustration of 
the 120 degree Neel order by next neighbor interactions. 
While some part of this spin liquid region in the phase 
diagram above kJc ^^^^ might be occupied eventually by 
magnetically ordered M-states, the existence of zero-flux 
spin liquid states for fairly large values of the quantum 
parameter is likely to persist, making the J1-J2 model an 
interesting candidate for spin liquid physics that deserves 
further attention. 



B. Effect of Ring Exchange 

In addition to the exchange interaction between a pair 
of spins, multi-spin interactions are also present in insu- 
lators with non-negligible charge fluctuations. While for 
spin one-half systems, the three spin interaction can be 
expressed in terms of two spin exchange, a new term is 
generated when considering four spin interactions. The 
four spin ring exchange term is: 

HBing = Jmng^{Pi+h.C.) (42) 

= h\sh<.h]j,iphlphkA^bu (43) 

where P4 is the four spin permutation operator, which is 
written out conveniently in terms of the Schwinger bo- 
son operators as above. The sites j, j, fc, I form a rhom- 
bus and a, P, "f, S d {"f, |} are spin indices. The above 
bosonic spinon representation of the ring-exchange may 
be checked by expanding it in terms of the spin opera- 
tors in the spin-1/2 case and verifying it has the required 
form'^*^. A more efficient way to prove this identity will be 
to establish a correspondence with the fermionic spinon 
representation of ring exchange which is done later in this 
section. 

From the t/U expansion of the Hubbard model^i, it is 
found that JRmg = 20^^/11^ and is positive. In the fol- 
lowing, when we consider ring exchange term we neglect 
the next-nearest-neighbor term, and vise versa. Evalu- 
ating the ring exchange term in the mean field ground 
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FIG. 8; (color online) The mean field phase diagram of the 
triangular lattice antiferromagnet with nearest neighbor Ji 
and next-neighbor J2 interactions, as a function of the ratio 
a — J2/J1 and quantum parameter k. The diagram is a 
quantitatively reliable solution of mean field equations in the 
region near to and below the red line (dashed and solid) where 
spinon condensation of the 7r-flux state occurs. Just above the 
solid red line, magnetic ordering at the Brillouin zone M-point 
occurs, while below it the vr-fiux spin liquid phase is obtained. 
The black line is the first-order phase boundary between the 
zero-flux and yr-flux spin liquid phases, while the blue line 
represents the onset of 120 degree magnetic order. The quick 
rise of this line to large values of k is due to frustration of 
the 120 degree Neel order by next neighbor interactions. If 
uninterrupted by magnetic order (of the M type) this may 
open a window of spin liquid for physical spin values in this 
model. 

state, one find that the result contains: 

{HRing)MF = SJRing{A^JA*k^klA*l, + C.c) + ... (44) 

which is the term we used to define fiux. The other terms 
are small, especially in the small-K regime. This term 
clearly has different signs for zero-flux and 7r-flux states, 
and favors the 7r-fiux states when Jumg is positive. Thus 
we also expect a possible transition from zero-flux state 
to TF-flux state when we tune ring exchange coupling. 

It is well known that the ring exchange term, de- 
rived from the Hubbard model, has positive (negative) 
coupling constant when expressed in terms of a per- 
mutation operator, if the loop length is even (odd)^. 
This can be seen as follows. Since the ring ex- 
change term arises from virtual electron hopping pro- 
cess lowering the energy, the effective term generated is 

-KL[c\„^CLcr,^c\^^_,CL-lar.-^ ' ' ' Clai + h-C] whcre 

CiCT are electron annihilation operators with site in- 
dex i and spin index ct, and > 0. Clearly this 
term includes a permutation operator for length L loop 

{-l)^-\-KL)c{^^Cia,cl^^C2a2 • • • C^^^_^ Ci^^ h.C. . So 
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for L even we have +Kl{Pl + h.c.). 

While the above discussion has been phrased in terms 
of fermionic spinons, we can now translate those observa- 
tions into the bosonic spinon language. It turns out that 
the permutation operator has the same form and sign in 
terms of both fermion (electron) operators and Schwinger 
boson operators. We give two arguments in support of 
this below. First, notice that spin 1/2 operators has the 
same representation in Schwinger boson and " Schwinger 
fermion" (electrons) scheme. Assume that spin indices 
a, f3 take values of -|-l(spin up) and -l(spin down). For 
S — 1/2 system we have cJ^Q^ = (25' + 2aS^)Sa,p + 
[(1 + a)S+ + (1 - a)S-] 5a^^i3 = bl^hp. Therefore we 
can simply replace fermion operator c by Schwinger bo- 
son operator b in the permutation operator. The sec- 
ond argument is to use Jordan- Wigner transformation 
to change electronic operators into boson. After replac- 
ing fermionic operators with bosons the Jordan- Wigner 
string operators becomes a constant under the constraint 
of constant (25) particles on each site. 

Now we evaluate the above ring exchange operator for 
longer length loops in the Schwinger boson mean field 
states. We focus on even length loops since we expect 
the dominant amplitudes to be given by the A.yS, which 
can only give rise to even length loops if they are the only 
amplitudes that count. The dominant term is 2 • 2^/^ • 
(yli,2v42,3 • ■ ■ ^L-i,L^L 1 + h.c). This is not exactly the 
term defining flux in the loop2^, which has an additional 
minus sign for each A*. Namely the result is 2 • 2-^/^ • 
(— l)^/^(term defining flux -I- h.c). It is then clear that 
for L = 4, 8, 12, . . . the ring exchange term favors 7r-flux 
in the loop, while for L = 6, 10, 14, . . . it favors zero flux. 
It is amusing to contrast this with the result for mean 
field theories based on a fermionic spin representation, 
where ring exchange favors states with zero-flux^^. The 
difference has to do with the different definitions for flux 
that naturally appear in these two theories. 

Small K Analysis: As argued previously, raising the 
value of JBing/Ji abovc a critical value, favors the tt 
flux state over the zero-flux state. We now study the 
effect of ring exchange term by evaluating it in the zero- 
flux and TT-flux states, in the limit of small k where an 
analytic solution may be obtained. Another benefit for 
going to small-K limit is that larger-ring exchange effect 
{e.g. length-six ring) is controlled by small parameter k 
as according to Tchernyshov, et. alm^, there- 

fore only the energetics of length-four ring, namely the 
rhombus, is important. 

The procedure described above is of course only the 
lowest order term, but we expect it to give us the correct 
qualitative behavior. For zero-flux state, the contribution 
from the ring exchange term is {k^ /6)JFiing per rhombus 
to lowest order. For 7r-flux state, it is — (k^/18) J^ing 
per rhombus to lowest order. Therefore, a positive Jmng 
clearly favors the 7r-fiux state. Since for every site, there 
are three bonds and three rhombi on average, we can find 



the transition point in the k — > limit to be : 

[Jmng/Jl]c = 1/8. (45) 

Several other theories predict a transition for the tri- 
angular lattice antiferromagnet while tuning this ratio 
jRing/Ji, although the phases appearing may be differ- 
ent. Motrunich's fermionic variational methodic predicts 
a transition from a " 7r-flux" fermionic mean field state 
to " zero-flux" mean fleld state when Jmng/Ji increases 
over about 0.35. A previous bosonic mean fleld phase 
diagramSS predict a transition from 120 degree Neel- 
ordered state to collinear-ordered state when this ratio 
increases over 1/3. It is interesting that we found a simi- 
lar transition between two disordered spin liquid state for 
sufficiently small n, with short ranged order at precisely 
these two wave vectors, as this ratio is increased. 

Numerical Solution: For finite k, the numerical solu- 
tion of the mean field equations have been obtained for 
the nearest-neighbor model below the critical value Kc 
for both states. We can make a lowest order pertur- 
bative treatment of next-nearest-neighbor and ring ex- 
change terms, namely we evaluate these two terms in the 
mean field ground state of nearest-neighbor model and 
thus obtain the transition point. We found approximate 
critical value of a = J2/J1 and JBing/Ji, listed in TA- 
BLE H] 



K 


Jmng/ Jl 


0.05 


0.126 


0.10 


0.127 


0.20 


0.128 


0.30 


0.129 



TABLE L Approximate finite k result for critical jRing/Ji at 
which the zero-flux state gives way to the yr-flux state on the 
triangular lattice. 



VI. SPIN LIQUID STATES ON THE KAGOME 
LATTICE 

The Kagome lattice has three sites in one unit cell, 
labeled by u, v, w (see FICEJ- 

We label a unit cell by r = (ri, r2) and label a site by 
(''I , f2)p where p = u,v,w. 

A. Projective Symmetry Group Analysis of the 
Kagome Lattice 

We would like to obtain the possible symmetric spin 
liquid states on the Kagome lattice described within the 
Schwinger boson approach. Again, a Projective Symme- 
try Group analysis will be carried out just as was done on 
the triangular lattice. A gauge transformation G is then 
described by three phase functions (f>u{^), 0i>(r), 4>w{^)- 
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FIG. 9: Coordinates system of Kagome lattice. Rhombus 
enclosed by dash lines is the unit cell at origin. Three sites 
with black dots form the basis of lattice. Black rhombus is 
an example of length-eight loops. 



For simplicity we will consider ansatz with Aij = 
^^aa' {biabja') Only. Thcsc amplitudes are naturally 
present for predominantly anti-ferromagnetic interac- 
tions, and can be strictly justified in large-N limit^^. 
Then, due to the frustrated nature of the Kagome lat- 
tice, the Invariant Gauge Group is Z2, generated by 

(pu{r) = (j)v{r) = (l)w{r) = tt 

thus we will study Z2 spin liquid states within the 
Schwinger boson approach. 

In the above coordinates system the space group is 
generated by two translations Ti, T2, reflection a and 
60 degree rotation . The first two generators preserve 
sublattices. The reflection exchanges u and v sublattices. 
The rotation is a cyclic permutation of three sublattices 
U ^ V, V ^ w, w ^ u. 

For each space group generator X we can associate a 
gauge transformation described by three phase functions 
'/>x,« (r) , 4>x,v (r) , (j)x,w (r) ■ 

The procedure of solving the for the different PSGs al- 
lowed by the algebraic constraints imposed by relations 
between symmetry elements is parallel to triangular case, 
except that one has to keep track of the index of sublat- 
tices. It turns out that the three phase functions can be 
chosen to be identical. The solution of algebraic PSG 
(details in Appendix Q is then identical to the triangu- 
lar lattice case except that the phase functions have one 
more index. 



<i>Tup{fi, ^2) 

'^T2,p(''l, f2) 
<^a,p(''l, ^2) 



(46a) 

Piirri (46b) 

P2n/2 + pinrir2 (46c) 

P37r +]3i7rr2(r2-l + 2ri) (46d) 



where p — u,v,w. Thus, in the case of the Kagome lattice 
too, there are not more than 8 symmetric Z2 spin liquid 
states, corresponding to the two values that each of the 
pi, p2, P3 can take on. 

If we now specialize to states that have non-vanishing 
nearest- neighbor amplitudes Aij^ we have more realiza- 



tions of the algebraic PSG. The only additional restric- 
tion(c./. SEC. mra is that bond (0, 0)„ ^ (0, 0)^ can be 
related to bond (0,0)u {0,0)w by both reflection and 
60 degree rotation. This gives one constraint p2 = l—pa- 
In all these ansatz the amplitudes Aij are real and of uni- 
form magnitude. 



B. Spin Liquid States on the Kagome 

The four realizations are listed below. We will use the 
same convention of k-space coordinates system as trian- 
gular case by taking the lattice constant as unity (the 
length of nearest- neighbor bond is 1/2). It is convenient 
to divide the four states into two pairs, two with pi = 
and two with pi = 1. The first condition implies that 
the generators of translation commute, the analog of the 
zero-flux state, while the second condition implies a flux 
of TT in a unit cell. It will turn out that the two zero-flux 
states correspond to the two states discussed earlier by 
Sachdev in the context of the Kagome lattice, while the 
two pi-flux states have not previously been discussed. 
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FIG. 10: (Color online) . k-space of Kagome lattice. The large 
hexagon is Brillouin zone of original Kagome lattice and the 
two pi = ansatz. The large rectangle is Brillouin zone of 
the two pi — 1 ansatz, which enclose a flux of vr in their unit 
cells. Red dots are the minima of the spinon dispersion for 
both pi = 1 states within the reduced Brillouin zone. Blue 
triangles are the locations of the minima of the lower edge of 
the two-spinon spectrum for both pi = 1 states. 

The Zero-Flux States: 

1. The first state we consider is characterised by 
Pi = 0, p2 = 0,p3 = 1. This is the large-N ground 
state identified by Sachdev for the nearest neighbor an- 
tiferromagnet, which he called Qi = —Q2 statei^. 

This state has zero-flux in both length-six hexagon 
and length-eight rhombus, hence we will also refer to 
it as [0IIex,0Rhom]. Its critical quantum parameter 
Kc ~ 0.54. The spinon dispersion has minima at cor- 
ners of Brillouin zone k = ±(27r/3, 27r/3), and so does 
the two-spinon spectrum. After spinon condensation it 
gives rise to the a/3 x \/3 magnetically ordered classical 
ground state. 
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2. The second state we consider is characterized by 
Pi = 0, P2 — 1,P3 — 0. This is another state considered 
by Sachdev and was called the Qi = Q2 statei^. 

This state has tt flux in length-six hexagon and zero- 
flux in length-eight rhombus, hence we will also refer 
to it as [TrHeXjORhom]. Its critical quantum parameter 
Kc ~ 0.50. The spinon dispersion has minima at cen- 
ters of Brillouin zone k = (0,0), so does the two-spinon 
spectrum. After spinon condensation it gives rise to the 
translational invariant classical ground state. 

The TT-Flux States: 

3. The third state we consider is characterized by pi = 
1, p2 = 0,P3 = 1. The ansatz is given in FIG.E|in the 
Appendix [HI and has not been considered before. Due to 
the presence of flux in the unit cell, a doubled the unit 
cell needs to be considered for the spinous. 

This state has tt flux in both length-six hexagon and 
length-eight rhombus, hence we will also refer to it as 
[7rHex,7rRhom]. Thus, in anearest neighbor model it is 
expected to have the highest energy of the four states for 
same k in the small k regime. 

Its critical quantum parameter is fairly large, Hc ~ 
0.93. The minima of spinon dispersion is at k = ±(7r/6 + 
nn, Tr/6 + mn), where n, m are integers(see FIG. [TUll . 
Thus the two-spinons spectrum lower edge has minima at 
±(7r/3 -|- nTT, 7^/3 + rmr) and (titt, rmr). In addition to the 
Brillouin zone center, corners and edge centers, minima 
are also present at the halfway points to the corners. 




FIG. 11: Ansatz of the pi = 1, p2 — 1, P3 = state. 
Presence of tt flux in the Kagome lattice unit cell leads to unit 
cell doubling, which is now defined by the dashed rhombus. 
Six sites with black dots form the basis. All nearest-neighbor 
Aij are real and of the same magnitude, and are positive in 
the directions shown. This spin liquid ansatz has zero flux in 
the length six hexagonal loops but n flux in the length eight 
loops such as the rhombus, hence [0Hex,7rRhom] . It has a 
large critical quantum parameter Kc = 2.0. 

4. The last state we consider is pi = 1, P2 = 1,P3 = 
0. The ansatz is given in FIG. and has not been 
considered before. 

This state has zero flux in length-six hexagon and tt 
flux in length-eight rhombus, hence we will also refer to 
it as [0Hex,7rRhom] . 

It is potentially the most interesting of the states 




kxa 
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FIG. 12: Two-spinons spectrum lower edge of the pi = 
1, J32 = 0, P3 = 1 state at K = Kc — 0.9. Axes are in di- 
mensionless unit kx,ya where a is lattice constant(two times 
nearest neighbor distance). Darker regions have lower energy. 



considered so far since its critical quantum parameter 
Kc ~ 2.0 is greater than unity. This means that this 
state is very likely to be a symmetric spin liquid state 
for spin-1/2 system if it can be realized. The large value 
of result of the dispersion of the relevant spinon 

bands being nearly flat along certain directions in the 
Brillouin zone. 

However, since this state has tt flux in the length-eight 
rhombus, it is, by the general argument of flux-expulsion 
by Tchernyshyov et. alm^ expected to have a higher mean 
field energy than Sachdev's ground state if we consider a 
pure nearest-neighbor Heisenberg model. In the small k 
limit the energy difference will be the order of k*. How- 
ever, it must also be kept in mind that mean field ener- 
getics is not precisely the same as the true energetics of 
the system. 

Given the interesting nature of this state we can ask 
what interactions are likely to stabilize it within mean 
field theory, relative to the Sachdev ground state. From 
our previous analysis of ring exchange we know that 
adding ring exchange interactions of the sign obtained 
from the Hubbard model favors zero flux in loops of 
length six and tt flux in loops of length eight. This is pre- 
cisely satisfied by loops in this ansatz as can be checked 
from the FIG. Therefore the Kagome antiferromag- 
net with ring exchange may be a good system to realize 
a spin liquid phase. 

The position of minima of spinon dispersion and two- 
spinons spectrum lower edge are the same as in the pre- 
vious state: pi = 1, P2 — 0,P3 = 1- However, the spinon 
dispersions of these two states are not identical. In par- 
ticular the present pi — 1, P2 = 1,P3 = state has 
three doubly-degenerate dispersion branches while previ- 
ously discussed the pi = 1, P2 = 0,p3 = 1 state has six 
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FIG. 13: Two-spinons spectrum lower edge of the pi — 
1, P2 ~ 1, Pa ~ state at k = Hc = 2.0. Axes are in di- 
mensionless unit kx,ya where a is lattice constant(two times 
nearest neighbor distance). Darker regions have lower energy. 
Although the two spinon minima occur at the same wavevec- 
tor locations as in the previous pi = 1, p2 = 0, pa — 1 state, 
the differences in the dispersion are apparent. 



non-degenerate dispersion branches. Other differences 
between these two states may be observed in the contour 
plots of the lower edge of the two spinon dispersions in 
FIG. 1131 and FIG. 1121 respectively. A variety of magneti- 
cally ordered states can arise from spinon condensation in 
both these cases; the magnetic states so obtained and the 
associated quantum transitions would be an interesting 
topic for future study. 



VII. CONCLUSIONS 

We have extended the projective symmetry group anal- 
ysis, previously used to classify fermionic mean field 
states, to bosonic mean field states of quantum antifer- 
romagnets on the triangular and Kagome lattices. This 
allowed us to access fundamentally different spin liquid 
states and ensure that potentially interesting states were 
not overlooked. On the triangular lattice a new Z2 sym- 
metric spin liquid state (tt-Aux state) emerged from this 
analysis and was studied. This state is distinct from 
the previously studied bosonic mean field state (zero- 
flux or Sachdev state), which can be seen in that they 
have different spinon and two-spinons dispersions . The 
nature of the ordered state resulting from spinon conden- 
sation turns out to be more intricate in this case, as is 
the quantum ordering transition, and is a topic for fur- 
ther study. Although this state is less favorable from the 
point of view of mean field energy as compared to the 
zero-flux state in the case of nearest neighbor interac- 
tions, it is stabilized by the inclusion of moderate next 



nearest-neighbor or ring exchange interactions with signs 
as derived from the Hubbard model. The mean field anal- 
ysis suggests that the triangular lattice model with next 
neighbor antiferromagnetic interactions might be favor- 
able for realizing spin liquid ground states. 

We have also studied Z2 symmetric spin liquid states 
on the Kagome lattice using the PSG technique. Again 
we find only a few candidate states, in particular two 
physically interesting states are found which are analo- 
gous to the TT-flux state on the triangular lattice. In par- 
ticular, one of these states has a remarkably large value 
for the critical quantum parameter Kc = 2.0 which im- 
plies that if realized, it would be stable against magnetic 
order for physical values of the spin. We find that ring 
exchange interactions would stabilize precisely this state, 
making the Kagome model with ring exchange an attrac- 
tive model to study in the search for spin liquid phases. 
Stability of these states to other forms of order, such as 
valence bond solid order, is important to understand but 
were not addressed in the present work. Also, the re- 
lation between fermionic mean field states, classified via 
the SU (2) PSGs and the bosonic spinon mean field the- 
ories studied in this paper with their U{1) PSGs remains 
unclear and is a topic for future study. 
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APPENDIX A: SOLUTION TO ALGEBRAIC 
PSGS OF TRIANGULAR LATTICE 

In this appendix we derive the allowed PSGs for sym- 
metry spin liquids on the triangular lattice character- 
ized by a Z2 gauge group. Essentially, we need to make 
sure that the constraints imposed by various identities 
between symmetry group elements are satisfied by the 
PSG. In particular the ansatz must satisfy the relations in 
equations l|16al to [T6li|) . It turns out that if a PSG satisfies 
all these constraints, then other relations between sym- 
metry group elements are automatically satisfied. This 
is because any group element written as a string contain- 
ing products of generators can be brought into a " normal 
ordered" form using just these operations where normal 
ordering implies a form of the kind: a'^ [iJiijT^^Tj^ where 
s e {0, 1}, r e {0, 1, . . .5} and h, ta £Z. This is suf- 
ficient to show no new constraints are imposed by other 
relations between symmetry operations. 

To simplify expressions we introduce forward difference 
operators Ai and A2, defined as Ai/(ri,r2) = /(ri -|- 
I,f2)-/(ri,r2) and A2/(ri,r2) = f{ri,r2 + l)-f{ri,r2). 

We should first consider how PSG changes if we do a 
gauge transformation G to the ansatz. After the gauge 
transformation, the ansatz will be invariant under G ■ 
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GxX ■ G-'^ = G ■ Gx ■ XG-^X-^ ■ X. Therefore, Gx 
should be replaced by GGx • XG'^X^^. In terms of 
phases, 

'/'x(r) ^ 0G(r) + </>x(r) - MX-\v)] (Al) 

As in Wen's derivation of fermionic PSG, we make 
the following assumption: using gauge freedom, we can 
make the ansatz explicitly invariant under translation Ti , 
or, 0Ti(r) = 0. For triangular lattice, this can be done 
by solving equations 0G(f'i,^2) + 4>Tt{ri,r2) - 4>Giri - 
l,r2) = 0. Since (l)x is a phase, unless explicitly men- 
tioned, all equations of phases in this section arc true 
modulo 27r. 

Then we add the generator T2. Since Ti-^TaTiTa"^ = 
X, following the procedure in the example in main text, 
we have the equation, 

Ai(?!)T2(ri,r2) = piTT 

where pi is a site-independent integer. Due to the 2tt 
periodicity, pi has only two distinct choices: and 1, 
and —pi = pi mod 2. Other integer parameters, p„ and 
p'„, in this section also have this property. 
Solution of (j)T.2 is then 

4>T2iri,r2) = (/)T2(0, ^2) +piTrri 

Here we make another assumption: using the gauge 
freedom, we can further make (^t'2(0i''2) = 0, while pre- 
serving (^Ti (r) =0. Then the solution simplifies to 

4>T2{ri,r2) = PiTrri (A2) 

In the above, we made two assumptions, c/jti (r) — 
and 0T2(O,?'2) = 0. These two equations can be satisfied 
if we have open boundary condition. However great care 
must be taken for periodic boundary condition before we 
proceed. In this section we always assume open boundary 
condition, and the two assumptions can be realized by 
exploiting gauge freedom. 

After this gauge fixing procedure we are still left with 
three gauge freedoms. The first one is adding a global 
constant phase, 

Gi : (/)i(r) = const. (A3) 

According to the equation ljAl|l this will not change any 
generator of PSG. We can use this freedom to fix one of 
Aij to be real positive. 

The second gauge freedom is 

G2: (I>2{ri,r2) =nri (A4) 

At first sight (j)Ti will be changed under this gauge trans- 
formation. According to HA1|I . after applying G2, 

0Ti(ri,r2) -> (t>2{ri,r2) + (l)Ti{ri,r2) ~ (l)2{ri - l,r2) 



However, we are free to add a site-independent constant tt 
to any 4>x{i^) because of the IGG structure. So this gauge 
transformation does not really change Gti and Gt2 ■ 
Similar to this one, we have the third gauge freedom, 

G3 ■■ 4iz{ri,r2) = 7rr2 (A5) 

G2 and Ga do not change Gti and Gt2 , but will cer- 
tainly modify other generators of PSG. Later we will use 
them to eliminate redundant parameters in our solution. 

We can now include generators of point group into con- 
sideration. First consider the reflection: cr. 

Algebraic constraints from Ti~^aT2(7^^ — I and 
r2"VTicr-i =1 are 

Ai(?!)o-(ri,r2) = P2^+pir27r 
A2(?!)o-(ri,r2) = p'^Ti+pirin 

Solution to these equations is 

(t>cr{ri,r2) = 0o-(O,O) + P2riTr + p'^r2TT + pirir2Tr 

Further constraint from aa = X is pj, = P3 mod 2 and 
20,(0,0) =P2^. 

Under the gauge transformation G2, the solution be- 
comes 

4'airi,r2) ^ (f>2iri,r2) + (f>airi,r2) - (f>2ir2;ri) 

<l>a{0,0) + iP3 + l){n+r2)7r+pirir2Tr 

while (^Ti and 0^2 do not change. Therefore we can al- 
ways assume that P2 ~ p'^ — mod 2. Finally we get 
the general solution l|17c|l in the main text. 

We are still left with two gauge freedoms at this point. 
Both of Gi and G2G3 do not change Gti, Gti, and Ga- 

Add the last generator to the system. Alge- 
braic constraints from Ti^^RilT2~^Ril^^ = 2 and 

^ 3 ^ 3 

T2~^ RilT2TiRil~^ =X are 

Ai0i?,, (ri,r2) = p'^Tr+pir2Tr 

3 

A2(?!'_R„ (ri,r2) = p'^n + pi{ri - r2 - l)Tr 
Solution to these equations is 

(r'i,f2) = (0,0) +P47rri -|-pg7rr2 

1 

+ Pirir2n + -pir2(r2 - Ijn 

Constraint from (i?-)^ = X is 12(/)h, (0,0) = 0. 

Further constraint from RiLaRiLU — X is p'/^ = p'^ = 
mod 2 and 4(/)fl'„ (0, 0) = 0. Therefore we can assume 

that (/>i?„ (0,0) =P37r/2- 

Under the gauge transformation G2G3, the solution of 
transforms to 

(?!>fix (?-i,r2) ^ yTT + (p5 + l)ri7r 

+ Pirir27r + \pir2{r2 - l)7r 
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while (?!iTi, 4>T2-: and (j)^ do not change. Therefore, we 
can always assume that P5 = mod 2. Then we get the 
general solution (|17d|l in the main text. 

A motivation for studying the range of states available 
comes from the numerical simulations^ on the triangular 
lattice S=l/2 model with ring exchange, where a state 
with a four fold degenerate ground state and no magnetic 
order was found in exact diagonalization studies. While 
this degeneracy is consistent with a gapped Z2 spin liq- 
uid, the quantum numbers of the degenerate states were 
not compatible with that expected for the simplest such 
spin liquid. The question then arises whether more com- 
plicated spin liquid states could reproduce the observa- 
tions. Although we answer this question in the negative 
for the TT-flux state, the nature of a phase that could pro- 
duce the observed degeneracies remains an interesting 
open question. 



APPENDIX B: PSGS OF ANISOTROPIC 
TRIANGULAR LATTICE 

Anisotropic triangular lattice can be realized in ex- 
perimental systems such as Cs2CuCl^. It is also pos- 
sible that the spin liquid state breaks triangular lattice 
space group symmetry down to anisotropic triangular lat- 
tice space group. A careful large-N study of Hubbard- 
Heisenberg model on anisotropic triangular lattice has 
been carried out by Chung, et. al^. They used a unit 
cell with four sites. Our PSG analysis can actually jus- 
tify their choice of unit cell within symmetric spin liquid 
states. 

We use the following coordinates system. 




FIG. 14: Coordinates system of anisotropic triangular lattice. 

This space group has also four generators. The first 
three are identical to the isotropic case: Ti, T2, and a. 
The last one is replaced by a 180 degree rotation i?^. 

: (r-i, ra) ^ (-ri, -ra) (Bl) 

The solution is identical to the isotropic case before 
solving the last generator. 

For the last generator i?^, from constaints 
Ti ^R^Ti ^ ^ = X and T2 ^R-kT2 ^R-k ^ — ^ ^ 



we have two equations 

Ai(/)i?,(ri,r2) =P57r 
A2(/'_R,(ri,r2) =P47r 

Solution to these equations is 

'f>R^{ri,r2) = 07?,, (0,0) + p^riTi + pir2TT 

Constraint from RttRtt = X is 40/?,, (0, 0) = 0, then we 
can assume that (0, 0) = P37r/2. 

Further constraint from aRT^aRj^ ~ 2 is p'^ = P4 
mod 2. 

However we cannot use gauge freedom G2G3 to elimi- 
nate Pi this time. 

The final solutions are 

0Ti(ri, ra) = (B2a) 

(l>T2i'ri, r2) = PiTTri (B2b) 

(t>criri, r2) = _P27r/2 + pirir27r (B2c) 

(f>R^i'ri, r2) = P3TT/2+p4Tr{ri+r2) (B2d) 

where pi, P2, P3, Pa are integers, either or 1. 

If nearest-neighbor A are nonzero, we have only two 
possible realizations. The difference from isotropic case 
is, this time not all A are real. Therefore we can, in 
principle, have time reversal breaking states. 

Two realizations: "zero-flux state": pi — 0, p2 — 
0, P3 1, P4 = 0. 

<l>T^{rl, r2) ^ (f>T2{ri, ^2) = (/)a(ri, ra) = 

(t>B^„{ri, r2) = 7r/2 

All B are real, diagonal A can be complex, A on the other 
two kinds of bonds are real (see FIG. I15|l . 




FIG. 15: Ansatz of zero-flux state on anisotropic triangular 
lattice. Diflterent arrows represent different amplitude of A. 
Single-arrow A are real. B are not represented in picture. 

"-TT-flux state": pi = 1, p2 — 1, P3 — 0, Pi — 1. 

<^Ti(ri, r2) = 

</'T2(?'i, r2) = Trri 

0<t(?'i, r2) = 7r/2 + rir27r 

0fl,(?'i, r2) = 7r(ri + r2) 

Diagonal A can be complex, diagonal B must be zero. A 
on the other two kinds of bonds are real, B on the other 
two kinds of bonds are pure imaginary (see FIG. 



19 




FIG. 16: Ansatz of vr-flux state on anisotropic triangular 
lattice. Different arrows represent different amplitude of A. 
Single-arrow A are real. B are not represented in picture. 



In both states, if diagonal A are neither real nor pure 
imaginary, we will have non-trivial flux (not zero or tt) 
through certain rhombi, then the ansatz breaks time re- 
versal symmetry. If we further impose time reversal sym- 
metry, diagonal A must be real in zero-flux state, and 
pure imaginary in 7r-flux state. However, if all B hap- 
pen to be zero, diagonal A can be either real or pure 
imaginary in any of these two states. 

The "tt-Aux" symmetric spin liquid state (with T- 
symmetry) on anisotropic triangular lattice, has the same 
spinon dispersion minima positions with isotropic case, 
therefore they will give rise to the same magnetic order- 
ing wavevectors (M-points: midpoints of B.Z. edges). 



APPENDIX C: PSGS OF KAGOME LATTICE 
1. Properties of the different mean field states 




FIG. 18: Ansatz of the pi — 0, P2 ~ 1, Pa = state which 
is the Qi = Q2 state in Sachdev's notation. On spinon con- 
densation this leads to the g = magnetically ordered state. 




FIG. 19: Ansatz of the pi = 1, P2 = 0, pa = 1 state. 
Rhombus enclosed by dash lines is the unit cell. Six sites 
with black dots form the basis. 



The zero-flux states: Ansatz for these states are shown 
in FIGs[n|and[Tl 




FIG. 17: Ansatz of the pi = 0, p2 = 0, pa = 1 state which is 
the Qi — —Q2 state of Sachdev which on spinon condensation 
leads to the \/3 x \/3 state. 



2. IT -Flux States For the pi — 1^ p2 = 0, — 1 state, 
the ansatz is given in FIG. 1191 

The mean field Hamiltonian for these two pi = 1 (tt- 
Aux) states is similar to that of 7r-flux state on triangular 
lattice, equations 1)27(1 and 1(29(1 . But the vector spinon 
field ^ has twelve components and the P(k) matrix be- 
comes a 6-by-6 hermitian matrix in the following form 



P(k) ^ JiAi 



U(k) C(k) 
,C(k) 6(k)^ 



where -4(k), /B(k), C(k) are 3-by-3 anti- hermitian ma- 
trices. 

The C(k) matrix has the same form for both pi = 1 
states. 



C(k) = 



_g-*fc3/2 
g-ifc2/2 




For pi — \, P2 = 1 , P3 = state the -4(k) matrix is 



/ 



^(k) = 



-ik3/2 



_gifc2/2 

2cos(/ci/2) 



_2cos(/ci/2) 

The B{k) matrix is 



/ 



6(k) 



_g-i'=3/2 __gifc2/2 

2isin(A:i/2) 



yg-ife2/2 2zsin(fci/2) 
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For pi = 1, P2 = 0, P3 = 1 state the ^(k) matrix is 



-1/03/2 



( ° 

e*'=3/2 2isin(fci/2) 

^_p-ik2/2 2isin(fci/2) 

The i3(k) matrix is 



S(k) 



/ gifc2/2 

-e*'=='/2 2cos(fci/2) 

^_^-tk2/2 _2cos(fci/2) 



The second and third gauge freedoms do not change 
Gti and ■ But they changes the other two generators. 
We win use them to fix the form of Gcr- The second one 
has no correspondence in triangular case, 

G2 ■■ 02,«(l") = +00 
02,t;(l") = -00 
02,tu(l") = 

where 00 is an arbitrary constant. The third one is sim- 
ilar to the gauge operation we used in triangular case. 



2. Solution to the Algebraic PSG on the Kagome 

The procedure for solving the algebraic PSG is de- 
scribed below. 

The algebraic form of the four space group generators 

is listed below. Note the unusual —1 in Ril. Generators 

3 

Ti and T2 preserve sublattices. 

Ti ■■ (r-i,r2)p ^ (ri -f l,r2)p 
T2 ■ {ri,r2)p iri,r2 + l)p 

a exchanges u and v, 

<y ■ {ri,r2)u {r2,ri)^ 
iri,r2)v {r2,ri)u 
{ri,r2)w -> {r2 

cyclicly permutates the three sublattices, 

-Rf : (r'i,r2)„ ^ (ri - r2 - l,ri)^ 
iri,r2)v {ri - r2,ri)^ 
iri,r2)w -> (ri - r2,ri)„ 

As before we can assume that 0Ti,p(r) = and 
0T2,p(0, 7*2) = 0, for p = u,v,w respectively. 

From Ti^^T2TiT2^^ = X, we have the equation 

Ai0T2,p(''i,'^2) =Pi7r 

for p = u,v,w respectively. Note that the constant pi 
must be the same for three sublattices. 
Solution is then 

0T2,p(''i,r2) =pi7rri 

At this stage we are left with more gauge freedom than 
the triangular case. The first one is a uniform rotation 
of boson phases on all sites. 

Gi : 0i,)9(i") = const. 

As before this guage transformation does not change any 
element of PSG. We will use this freedom to fix one of 
Aij{e.g. bond (0,0)„ — (0, 0)tu) to be real positive. 



G, 



where p = u,v, w. 

The fourth gauge freedom does not change Gti, Gt2 

and Go-. But it changes Gr^. We will use this freedom 

3 

to fix the form of Gjj^ . 

3 

G4 : 04,«(?'i, r2) = 7r(ri + r2) 
(t)4:,v{ri,r2) = 7r(ri + r2) 
(f>4.,w{ri,r2) = 7r(ri + r2 + 1) 

We can now introduce point group generator a. 
Algebraic constraints from Ti~^aT2a~^ = X and 
Ts^VTict^i =X are 

Al0o-,p(''l,''2) = pi^'K+P\T2'K 
A20a,p(fl,f2) = PaTT+piriTT 

where p = u,v,w, p'^ and p'-^ are integer constant inde- 
pendent of unit cell index (ri, r2) and sublattice index p. 
Solution to these equations is 

0<T,p(?'i,?'2) = 0(T,p(O,O) -f P2^i7r + p3r27r-|-pirir27r 

However, we have no a priori reason to say 0cr,p(O, 0) are 
independent of sublattices index p. 

Further constraint from gg — X in p'2 — p'^ mod 2 and 

20.,I„(O,O) = P2TT 

0^,1.(0, O)-f0^,„ (0,0) = P27r 

where p2 is an integer constant. This fixes 0cr.u,(O, 0) to be 
P27r/2 but leaves one freedom for 0cr,„(0, 0) and 0cr,i,(O, 0). 
As in triangular case we can use G3 to make = P3 = 0- 
Because after applying G3 the solution to (f>„ becomes, 
according to HAip . 

0<T,p(?'i,?'2) 0CT,p(O,O) + {p'2 + l){ri + r2)7r +pirir27r 

We now apply G2. 

0<T,«(7-i,r2) 0CT,„(O, 0) + 200 +pirir27r 
0cr,i,(?'i,r2) 0<j,t,(O, 0) - 200 +pirir27r 

0a,«,(?'l,?'2) 0<T,«;(O,O) +pirir27r 
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Therefore we can always make (j)cr,u = 4'cr,v = P27''/2 by 
choosing appropiate 0o- Then the three phase functions 
(/)ct,p are identicaL 

Add the last generator to the system. Alge- 
braic constraints from Ti~^ RilT2~^ R^~^ = I and 
T2^^ R2lT2TiRil^^ =X arc 

A2(l)R„,p{ri,r2) = PgTT +pi(ri - r2 - l)7r 
Solution to these equations is 

^fl^,p(?'i,?'2) = (t)R^.p{^,0) +Pi'^ri +p'r,Trr2 

1 

+ Pirir2iT + -pir2[r2 - IjTr 
Again we should not assume that (f>R^ p(0, 0) are inde- 

3 

pendent of sublatticcs index p. 

Further constraint from R^aRiia = I is p'^ ^ p'^ = 
mod 2 and 

20fl^,,(O,O) = 

3 



Therefore we have i,(0, 0) = P37r/2. 

A similar constraint from aRi^aRi^ = I is p'^ = p'^ = Q 
mod 2 and 

2(/'fl,,„(0,0) = P47r 
'/'fli ,i>(0, 0) + = p47r 

3 3 ■ 

Combining with previous constraints on (f)R^ „ we con- 

3 

elude that p4 = py, mod 2 and 4)r^ = p37r/2. 

Under the gauge transformation G4, the solution of 
(l)R^ transforms to 

0i?^,p(?'i,?^2) ^ y7r + 7r + (p5 + l)ri7r 

+ Pirir27r + ^Pir2{r2 - 1)tt 

while (pTi.pi 0T2,pj and (/)cr,p do not change. Therefore, 
we can always assume that = mod 2. The result- 
ing constant tt can be neglected because it is an IGG 
operation. Then we get the general solution (|46d|l . 
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